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Abstract 

Electrical properties of metal-semiconductor (M/SC) and metal/oxide/SC structures built with 
Zr or Zr02 deposited on oxygen-terminated surfaces of (OOl)-oriented diamond films, comprised of 
a stack of lightly p-doped diamond on a heavily doped layer itself homoepitaxially grown on a Ib 
substrate, are investigated experimentally and compared to different models. In Schottky barrier 
diodes, the interfacial oxide layer evidenced by high resolution transmission electron microscopy and 
electron energy losses spectroscopy before and after annealing, and barrier height inhomogeneities 
accounts for the measured electrical characteristics until flat bands are reached, in accordance with 
a model which generalizes that of R.T. Tung [Phys. Rev. B 45, 13509 (1992)] and permits to 
extract physically meaningful parameters of the three kinds of interface: (a) unannealed ones, (b) 
annealed at 350°C, (c) annealed at 450°C with characteristic barrier heights of 2.2-2.5 V in case (a) 
while as low as 0.96 V in case (c). Possible models of potential barriers for several metals deposited 
on well defined oxygen-terminated diamond surfaces are discussed and compared to experimental 
data. It is concluded that interface dipoles of several kinds present at these compound interfaces 
and their chemical evolution due to annealing are the suitable ingredients able to account for the 
Mott-Schottky behavior when the effect of the metal work function is ignored, and to justify the 
reverted slope observed regarding metal work function, in contrast to the trend always reported 
for all other metal-semiconductor interfaces. 
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I. INTRODUCTION 


Metal-semiconductor interfaces are necessary to implement sensors and electronic devices, 
and diamond does not escape to this need. Building well defined interfaces at nanometer 
scale is specially necessary for electrical rectifiers, which can gain high benehts from the bulk 
properties of boron doped diamond like its very high thermal conductivity, hole mobility and 
electrical breakdown held [l, 2| if specihc properties of the interface are also fulhlled. These 
are mainly good adhesion, thermal and chemical stability at elevated temperatures, com¬ 
patible with the superior possibilities of diamond, and potential barriers ensuring both low 
electrical losses under forward voltage and minimal reverse currents even at high tempera¬ 
tures. Several attempts have been proposed, relying either on materials which react with the 
Dare surface of diamond to form carbides, like silicon or refractory metals and compounds 


^-l6|, even hetero-epitaxially grown on diamond [^, or metals on either the oxygenated or 
hydrogenated surface of diamond j8-12|. On bare diamond surfaces, potential barriers close 
to 2 V, stable at temperatures higher than 500°C have been obtained but reverse current 
densities and ideality factors were much greater than expected from the thermionic mech¬ 
anism alone at least at room and moderate temperatures, and some of these contacts were 
ohmic. These poor rectifying current-voltage characteristics have been assigned to inho- 


iQ- 


mogeneous and defective interfaces, specially with carbide forming metals [8|, IIJ . Lower 
potential barriers are implemented on the hydrogen-terminated diamond surfaces due to the 
dipole H+'^—0“"^ responsible for the negative electron affinity of these surfaces, where 'd is 
the averaged fraction of the elementary charge per atom. But the chemical stability is not 
guaranteed at temperatures as high as in other interfaces and reverse current densities turn 
out not to be lower than 10“^ A/cm^ at room temperature, still higher than the thermionic 
limit. Improvement has been achieved with the help of an oxide layer deposited on the H- 
terminated surface to build held effect transistors [ij, ll5| but this solution is not relevant for 
Schottky barrier diodes hrstly because it is very difficult to control accurately and in a repro¬ 
ducible way an oxide interlayer made of foreign species with thickness in the subnanometer 
range. Secondly, the Schottky barrier heights are known to be much smaller than 1 V and 
even zero for noble metals on the hydrogenated diamond surface, thus inducing far too 
high reverse currents in rectihers. Because of an inverted electric dipole at the oxygenated 
surface of diamond, the largest potential barriers are generally obtained with either noble 


3 















(Au) or transition metals deposited on this type of surface [9j, but some authors have shown 
that thermal treatments up to 500 or 600°C are able to decrease the barrier height down to 
1.2 V [^, 18|, probably because of the cancellation of the electric dipole still 

preserving a good adhesion of the metallic layer. More recently, we have shown that a less 
electronegative and more easily oxidizable metal like Zr was able to reach even more beneh- 
cial properties for the Schottky contact, namely an even lower potential barrier height while 
maintaining the reverse current density after an anneal at 450 °C close to 3 x 10“^° A/cm^ 


at room temperature 


19| and an ideal current-voltage behavior for most diodes. However, 


a comprehensive picture of metal oxygenated-diamond interfaces did not emerge clearly in 
the last hfteen years, neither for the transport properties of carriers through interfaces, nor 
for the question of how the potential barrier height is determined, contrary to the case of 
hydrogen-terminated diamond surfaces Ig. In the present paper, mechanisms of carrier 


flow and signihcance of barrier heights are investigated on the one hand with more details 
and analysed with the help of several models. 

On the other hand, the set of works quoted above demonstrates that both chemistry and 
physics must be both considered when dealing with Schottky barriers. This is specially true 
in the case of diamond where a large dispersion of barrier heights occurred for any metal 
js] in literature. Much work was devoted in the past to general models able to explain and 
predict band misalignments at interfaces compounded of various materials, and a recent 
review of the various approaches, including a mixed physical-chemical point of view, has 


been published by R. T. Tung 


2n| . But, as often the case for more than one decade in 


the literature excepted in ref. Iq, no experimental data regarding diamond interfaces has 


been brought and discussed. The present study aims at hlling this gap, giving evidence of 
the presence of an oxide interfacial layer and discussing the relevance of various models of 
barrier height regarding data presented in this article and others available in literature for 
oxygenated diamond interfaces. In the second section, a detailed study by high resolution 
transmission electron microscopy (HRTEM) and electron energy loss spectroscopy (EELS) 
is developed to show the characteristics of the interfacial layer at a sub-nanometer scale, 
the more noticeable being the systematic presence of oxygen at interface before and after 
annealing. In the third section, a hrst level analysis of the current-voltage characteristics 
relevant to the depleted semiconductor interface is presented and representative parameters 
of the junction are derived, conhrming the presence of an interfacial layer in-between Zr 
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and diamond and the large decrease of the barrier height after annealing. In the fonrth 
section, potential barrier inhomogeneities at different scales are analyzed and evidence of 
positively charged centers in the interfacial oxide is given. The conseqnences npon cnrrent- 
voltage characteristics are discnssed and it is demonstrated that the jnnctions experience a 
change from depletion to accnmnlation regime at snfficient forward bias. Two models able 
to take into acconnt the inflnence of the barrier height inhomogeneities are developed and 
matched to experimental data, allowing to derive meaningfnl parameters of each interface 
which indicate that large barrier inhomogeneities prevail before annealing while they almost 
completely vanish after snfficient annealing. The previons resnlts are nsed in the hfth sec¬ 
tion to address the problem of whether snch an accnrate control of the interface and metal 
choice can allow to tailor and predict the potential barrier heights of Schottky jnnctions 
on diamond. Gathering all the data available for well characterized interfaces of varions 
metals on oxygen-terminated diamond permits to nnveil a completely nnnsnal relationship 
between potential barrier heights and metal work fnnctions. Modeling this new resnlt needs 
to consider the contribntion of all electrical dipoles which may be present at these componnd 
interfaces. Conclusions are drawn in the sixth section. 


II. NANOSTRUCTURE OF ZIRCONIUM ON OXYGEN-TERMINATED DIA¬ 
MOND INTERFACES 


Schottky junctions are implemented by electron-gun deposition in ultra high vacuum of 
Zr on a stack comprising a lightly boron doped homoepitaxial diamond layer grown on a 
heavily doped one, itself grown on a Ib substrate. The surface of the lightly doped layer 
is submitted to a photo-chemical UV-ozone treatment for two hours at room temperature, 
so that oxygen terminations prevail on the diamond surface. Other details of the diamond 
growth and interface preparation are given in reference 19|]. Three types of sample are 
elaborated, the hrst one being annealed at temperatures not exceeding 300°C, while the 
second and third ones are annealed respectively at 350°C and 450°C; and the respective 
interfaces are labeled (a), (b) and (c). Identihcation of the oxide layer and evaluation of 
its thickness are carried out by a combination of high resolution transmission electronic 
microscopy (HRTEM), electronic energy loss spectroscopy (EELS), and high angle angular 
dark held-TEM (HAADF-TEM) techniques. HRTEM, HAADF and EELS measurements 


5 



were performed using the beam of a JEOL 2010F STEM microscope directed on electron- 
transparent specimens of the interfaces labeled (a) and (c), fabricated by a lift-off technique 
in a focused ion beam (FIB) FEI Quanta 200-3D dual beam microscope. 

The characterization of the oxide layer have been carried out in three steps: chemical 
contrast TEM (CTEM) identihcation of the metal/oxide/diamond interface, scanning TEM- 
EELS study of the previously identihed oxide layer and HRTEM imaging of the interface. 
In a hrst approach, bright-held (BE) images in CTEM mode are used to identify the inter¬ 
face. BE imaging variations in the electron absorption and dependence on the diffraction 
angles allow to distinguish between diamond and Zr material. In a second step, EELS 
spectra are acquired to chemically identify the oxide layer. During EELS spectra acquisi¬ 
tion, HAADF-TEM operation mode is used to identify the diamond/Zr interface, insofar 
contrasts in HAADF imaging are due to different Z-number of the elements in the sample. 
Secondly, EELS prohles are carried out, revealing the oxide layer through the observation 
of the EELS oxygen-related peak (see Fig. [TId). Thirdly, once CTEM and EELS techniques 
have been applied, HRTEM images of the previously identihed oxide layer are acquired. 
Fig. [1^ shows EELS spectra recorded across the diamond/Zr interface after the thermal 
treatment, at the positions indicated in the inset. Presence of oxygen is revealed at the 
diamond/Zr interface (position C). Moreover, the simultaneous observation of the Zr-M2 to 
4 and the 0-K peaks at position C discloses the formation of an oxide layer at interface, 
however slightly diherent from bulk Zr02 since the peak Zr-M5 is missing. In one side of 
the interface, a diamond-related shape of the carbon core-loss is evidenced, while typical Zr- 
related peaks are observed in the other side. Fig.[T]D shows normalized intensity prohles of the 
oxygen-related peak (see spectra at position C) acquired across the diamond/oxide/metal 
interfaces before and after thermal treatment. To improve signal/noise ratio, ten spectra 
were averaged. The latter prohle exposes how oxygen distribution changes after thermal 
treatment, showing a real oxide thickness reduction (FWHM varies from 0.6nm to 0.3nm) 
while the amplitude is increased by a factor 1.5, suggesting an augmented oxygen density. 
Note that the FWHM of the EELS linescan results from the convolution with the probe size 
of 0.2nm. Simultaneously, the EELS analysis of the oxygen content in the Zr layer indicates 
an averaged concentration of 5 % in the as-deposited metal while it decreases to 1 % after 
annealing. Finally, analysis of HRTEM images like shown in the inset of Fig. (Tb indicates 
some in-depth roughness in the interfaces of kind (a) whereas a smoother interface takes 
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Depth (nm) 


FIG. 1. a) EELS spectra recorded in an interface of kind (c) along the linescan drawn in the 

HRTEM image shown in the inset. The peak Zr-M2 is also visible near 345 eV but Zr-M5 normally 

detected in a 10 nm thick Zr02 layer near 190 eV is missing in the interfacial layer; b) Normalized 

intensity profile of the oxygen related peak across the Zr/diamond interface before (grey curve) 

and after annealing at 450°C for 30 minutes (black curve), corresponding to interfaces of kind (a) 

and (c) as labeled respectively in the text. In the inset, HRTEM filtered images of the Zr/oxygen- 

terminated diamond interfaces of kind (a, unannealed) and (c, annealed at 450°C) are shown. 

Superimposed contrasts of diamond and Zr can be observed in the Zr/diamond interface of kind 

(a), with white arrows to guide the eyes through the diamond/Zr interface. 

7 

























place after annealing in contacts of kind (c), leaving however an abrnpt transition across the 
interface between carbon atoms and those of the interlayer in both cases. 

Titaninm and Zirconium are well known to incorporate oxygen and possibly nitrogen when 
they condense onto a substrate after evaporation in vacuum, an effect which is clearly evi¬ 
denced by a pressure decrease by about one order of magnitude inside the deposition cham¬ 
ber. Consequently, it is not surprising to find 5 % of oxygen in the as-deposited Zr layer. 
Nitrogen is not detected. Increasing the temperature leads to augmented internal pressure 
and elastic energy of the Zr lattice, so that it becomes less expensive in term of total internal 
energy to form an oxide at the metal interfaces, both because the internal pressure is de¬ 
creased and the enthalpy of oxide formation is negative. Such effects explain why the oxygen 
content of the Zr layer is decreased after annealing. More paradoxical is the sharpening of 
the oxide interlayer shown by EELS spectra because an increase of the oxygen content at 
interface should be expected from the evolutions described previously. At least three reasons 
can reconcile this apparent contradiction: (i) the beam spot (with an approximate size of 
0.2nm) averages in-depth roughness during the EELS capture and, as a consequence, the 
recorded 0-K profile of Fig. [lb shows a widened profile in the interface of kind (a) before 
thermal treatment; (ii) we know that oxygen termination coverage has a lower density than 
the last carbon atom layer of the diamond lattice even after the most efficient treatment 


like UV-ozone oxidation used in this study 


2l|, a fact which correlates very well with the 


previous statement since a natural non-uniformity results before annealing; (iii) as demon¬ 
strated in the following sections, a strong change in electronic properties of the interface 
follows the 450°C anneal, correlated with the electron affinity variation of the diamond sur¬ 
face induced by the cancellation of the oxygen terminations. However, because oxygen is still 
present at interface after annealing, one is forced to conclude that bonding has been deeply 
modified, with a probable transformation from a physisorbed metal layer on the oxygen 
terminated diamond surface likely through the hydroxyl terminations (see subsection II V All . 
towards stronger bonds between carbon and oxygen atoms, characterized both by shorter 
bond lengths and a higher and more uniform density. Taking into account the convolution 
with the probe size, the oxygen thickness detected by EELS after 450°C anneal in interfaces 
of kind (c) suggests that the oxide interlayer comprises two planes of oxygen atoms sepa¬ 
rated by a Zr atom plane, since the distance between two oxygen planes is 0.256 nm in cubic 
zirconia, thus defining a full unit cell of Zr02 as the whole interlayer thickness. Next sections 





are devoted to determine and understand how oxygen species at the interface influence the 
electrical properties. 


III. CURRENT-VOLTAGE CHARACTERISTICS OF THE ZIRCONIUM-OXYGEN- 
TERMINATED DIAMOND INTERFACES 


A. Theoretical background of the relationship between current and voltage 


Current-voltage characteristics and mechanisms of carrier transport through the potential 
barrier present in metal-semiconductor junctions have been discussed extensively in review 


articles and textbooks 


2Q|, 


22|-|24| ■ The z direction being perpendicular to the x — y plane of 


the interface, the current can be calculated in the semi-classical framework of the Maxwell- 
Boltzmann distribution of majority carrier velocities, by : 


h. = 



Co Vz{E) D{E) dE dx dy 


( 1 ) 


Co being the elementary charge, Vz{E) the carrier velocity in the z direction and D{E) the 
majority carrier distribution just at energies E above the point where the potential barrier 
culminates. If the barrier is homogeneous, at a constant potential throughout the x-y 
plane, the current density can be written from the previous expression, in the case of a 
voltage Vj applied to the junction [22|, l25| : 

eoNjiiVcoii 




■exp {-eo^B/ksT) [exp(eoUj//csT) - 1] 


( 2 ) 


1 T Xcoll /Vdif f 

with Nm being the effective density of states at the majority carrier band edge, ks the 
Boltzmann constant, T the absolute temperature, Vcoii the collection velocity at interface 
due to the thermionic effect and equal to one quarter of the thermal velocity, m* 


2nm* 


the conduction effective mass of majority carriers, and Vdiff their diffusion velocity. The 
previous expression takes into account the two transport mechanisms of majority carriers 
which work in series, due respectively to the thermoionic velocity at interface and diffusion 
velocity inside the depletion zone, whereas it neglects bulk recombination. The diffusion 
velocity simplifies into the product of the majority carrier mobility by the electric field at 
interface if the band bending exceeds some ksT/eo units and is generally not the limiting 
factor for the direct current, as conhrmed in the samples under study from the comparison 
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of the two velocities. The numerator of the hrst factor in Eq. ([2]) can be rewritten A*T‘^ 
where A* = (Aneok'^mo/h^) x {m*/mo), the hrst quantity inside the parenthesis being the 
Richardson constant for emission of electrons into vacuum with a mass mg at rest, Aji = 
120 A cm“^ K~^. Because Vdi/f decreases with the increase of the forward bias voltage 
Vj, the pre-exponential factor of Eq. ([2]) becomes voltage dependent. Other reasons for 
such an effect involving also the factor exp (—eo$B/A;sT), are (i) the image force lowering 
of the barrier, (ii) tunneling through the top of the barrier possibly assisted by the electric 
held, (hi) voltage drop inside an interfacial layer between the metal and semiconductor, 
(iv) charging interface states in equilibrium with the semiconductor and (v) a spatially 
inhomogeneous barrier as discussed further. However, the second ehect is noticeable only 
in heavily doped semiconductors, a case which is not achieved in this study since the boron 
concentration of the lightly doped layer is close to Na = 5 x lO^^cm”^. As one almost 
always wants to characterize the Schottky contact with a single potential barrier, the more 
common approximation consists in keeping only the zero bias voltage barrier height 
in the hrst exponential factor of Eq. ([2]) and gathering all the voltage dependent ehects 
into the last factor. Since all these ehects contribute to a slope of the logarithm of the 
current which is smaller than e^/ksT, an ideality factor n larger than one is introduced 
in the denominator of the exponent of exp(eoV^/fcBT) and often considered as a constant 
at least for a given temperature. Such an hypothesis is strictly valid in the case of an 
interfacial layer where the voltage drop is a constant fraction 1 — 1/n of the band bending 
U = ~ Vj ~ Vdop (with coVdop being the energy diherence between the Fermi level and 

majority carrier band edge in the bulk semiconductor) mainly determined by the oxide 
capacitance and density of interface states 2^ when mechanisms (ii) and (iii) are combined. 
But it is not adequate if the interfacial layer either allows tunneling of carriers as the main 
transport mechanism or does not display any additional barrier for majority carriers because 
of band alignments similar to those occurring in type II heterostructures, or in the case of 
barrier inhomogeneities and even with image force lowering. Other effects, which are due to 
quantum reflexion, transmission and phonon scattering of carriers at the interface are voltage 
independent and generally very weakly temperature dependent so that it is convenient to 
lump them into an effective Richardson constant A** which is weaker than A* by some tens 


of percent in the case of intimate Schottky contacts 


3 . 


yielding the standard value of 


90 A cm ^ K ^ usually considered for diamond. This effective Richardson constant can be 
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further severely decreased because of the attenuation by a factor a of the collection velocity 
Vcoii due to additional scattering of carriers by charged centers and roughness, and transport 
mechanisms through the interfacial layer when they exist. Finally, the non ideal current 
density is rewritten as : 




aA**T‘^ exp 


V ksT J 



( 3 ) 


Experimental current-volt age characteristics can then be accounted for by such an ex¬ 
pression with the ideality factor n as an adjustable parameter, provided aA** is known and 
the series resistance Rg is adjusted in order to £t the additional voltage drop RgR appearing 
in the total measured voltage V = Vj + RgR- This is done in Fig. [2] for the three interfaces 
near room temperature. Two kinds of deviation between experimental data and the theo¬ 
retical laws may appear as it can be seen in Fig. |2] ; the hrst one is due to a non constant 
ideality factor as long as the voltage drop in the series resistance is negligible (typically below 
0.5 A/cm^ in these diodes at room temperature) and the second one is due to less common 
transport mechanisms which are discussed in the following. It is worthwhile to notice that 
the rectification ratios are more than twelve decades and that no generation-recombination 
current is detected on these current-voltage characteristics, which would have involved a 
symmetrical behavior of the absolute value of the current around the origin at low voltages. 
Such an almost ideal behavior strongly suggests that the oxygen terminated surface of dia - 
mond passivates very efhciently the interface defects able to develop leakage currents 


The determination of aA** needs temperature variations in such a way that the saturation 


current density over T^, namely a A** exp 


fcsT 


, obtained from the extrapolated value 


of the current-volt age characteristic at Vj = 0, can be drawn. Such Richardson plots are 
displayed in Fig. [3l The effective Richardson constants aA**, ideality factors at room tem¬ 
perature and parameters described in the following are indicated in the Table HI It must be 
stressed firstly that the aA** values are at least ten times lower than the standard value 
A** and secondly that a careful inspection of the Fig. [3] shows that the data points follow 
systematically a slight upward curvature, indicating a barrier height increase with temper¬ 
ature. The hrst result conhrms indubitably the presence of an insulating interfacial layer, 
which is able to increase considerably the quantum rehection coefhcient and decrease the 
effective velocity of carriers, mainly because of scattering due to roughness, chemical disor¬ 
der and differential dielectric response. The second one is the preliminary indication of the 
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Applied voltage (V) 


FIG. 2. Current density as a function of applied voltage at 300 K in the three interfaces made of 
Zr deposited on oxygen terminated diamond (a) with annealing temperature not exceeding 300 °C; 
(b) annealed at 350 °C; (c) annealed at 450 °C. Dash-dot lines are calculated with the help of Eq. 
(f3|) , a A** values included in table |I] which are obtained from fits achieved in the two following 
figures, and adjusted constant values of the apparent barrier height and ideality factor in each 
case, that are respectively: (a) 1.83 V and 1.57 ; (b) 1.47 V and 1.29 ; (c) 0.97 V and 1.05. Voltage 
takes into account the voltage drop in a constant series resistance of 600 D in case (a) and 300 D 
in cases (b) and (c), in good agreement with the theoretical ones, but current densities above a 
few A/cm^ turn out to be higher because of the conductivity increase due to hole injection from 
the heavily doped layer. Full lines rely on models which incorporate barrier height lowering and 
are described in section HVl with Eq. ([7]) for the cases (b) and (c), while Eq. (flU) is used for the 
case (a). Horizontal arrows mark the flat band situation in each junction. The effect of image 
force lowering is displayed in the inset for the case (a). In the hatched area, self-heating becomes 
important. 
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FIG. 3. Richardson plot for the three same interfaces as in Fig. [2l Averaged barrier heights 
deduced from the slope of each straight line are respectively: (a) 1.84 V; (b) 1.47 V; (c) 0.95 V; 
very close to that used in Fig. [2l 


presence of barrier inhomogeneities, a conclnsion similar to that obtained in ref. 29| abont 
silicides-silicon interfaces, also implying additional conseqnences which have to be analyzed 
more thoronghly. 
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Sample type and 

aA** 

n 

B,oo 

0-$ 

annealing temperature (°C) 

A cm“^ K~^ 

(at 300 K) 

V 

mV 

(a) < 300 °C 

11 

1.52 

1.93 

71 

(b) 350 °C 

2.8 

1.29 

1.57 

87 

(c) 450 °C 

8.5 

1.00 to 1.07 

0.93 

54 


TABLE I. Parameters of the Schottky barrier in the three interfaces according to the model of 
section imi 


B. Richardson constant, ideality factor and potential barrier behaviors 


Zero voltage barrier heights and averaged ideality factors are displayed in Fig. 0] as a 
function of temperature. On the first hand, the former quantities experience indeed an 
increase with temperature with generally an asymptotic behavior towards highest temper¬ 
atures, while the latter ones are decreasing, starting from very different values at room 
temperature (which are reported in Table [T|) and displaying very different negative slopes. 
As noticed previously, such variable ideality factors cannot be ascribed to a voltage drop 
within an interfacial layer bearing an additional barrier to that built in the semiconductor, 
because they would not be coherent with oxide interfacial layers which still exist at the three 
interfaces. On the other hand, experimental data plotted in Fig. [2] deviate from the simu¬ 
lated curves drawn in dash-dot lines on the basis of a constant ideality factor when voltage 
drop in the series resistance and other effects remain negligible. A hrst attempt to £t more 
accurately experimental results uses an alternative dehnition of the current density, which is 
also more physically meaningful. It consists in using a temperature and voltage dependent 
barrier height <Fb(T, V^), so that the current density can be rewritten more usefully for the 
homogeneous barrier case: 


>^.11 


aA**T^ exp 


/ eo$^(T,V,) \ 

V ksT ) 



(4) 


Introducing an inhomogeneous barrier which follows a Gaussian spatial distribution in the 
X — y plane of the interface, characterized by an average value <I>b(T, V^) and a standard 
deviation cr$(V^), permits to derive 


M 


29| a law similar either to Eq. (jS]) or (0]), with an 
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effective barrier given in the second case by : 


^BiT,V,) = ^BiT,V,) 


2kBT 


( 5 ) 


For each temperature, it is always possible to convert Eq. (|1]) into Eq. (|3]) through the 
dehnition of a voltage independent barrier height, conveniently taken at = 0, and a 
voltage dependent ideality factor n{T, Vj), valid in all situations, the inverse of which being 
evaluated from the derivative of the logarithm of the current density given by Eq. ([3]) and 
(jlj) when the reverse saturation current is neglected: 


(i(ln Jz 


^ ^ _ d^B{T,Vj) ^ ^ ^ d^B{T,U) 


( 6 ) 


eo dVj n(T,Vj) * dVj * ' dU 
However, in this hrst level model, no explicit dependence of the parameters is given as a 


function of junction voltage Vj or band bending U. Some cases of potential barrier distribu¬ 
tions have been analyzed analytically by R. T. Tung 30|, who showed that the dependence 
of the parameters of Eq. ([5]) was not linear with the band bending and hence with the 
junction voltage, on the basis of a physical representation of the barrier lowering by a £c- 
titious dipole located at interface. An additional source of non ideality may be present 
if a spreading resistance effect occurs, involving inhomogeneous current density inside the 
semiconductor because of current crowding near the pinch-off patches with a lowered barrier 
and located close to the interface j^, 32|. This effect is normally expected to be relevant at 
current levels where the series resistance voltage becomes larger than some fcsT/eo. But the 
combination of multiple barriers and incomplete ionization of doping impurities may induce 
deviation from the ideal current-voltage characteristic at much smaller currents because of 
severe weakening of the carrier concentration and narrowing of the current paths in the low¬ 
est barrier patches which cannot be assessed easily from the comparison to a voltage drop in 
an averaged series resistance. This kind of non ideality will be discussed in the next section. 

In the simple model presented above, the barrier fluctuations are characterized by their 
average and standard deviation but are not supported by any physical model, thus preventing 
the derivation of any relationship between physical properties of the interface and electrical 
ones. Moreover, the behaviors of the barrier lowering and ideality factor with voltage and 
temperature cannot be deduced. Therefore, the reasons for these barrier inhomogeneities 
have to be analyzed more accurately and an analytical model which will generalize that 


developed by R. Tung 


30| will be implemented in the next section, concurrently with a 
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Temperature (K) 


Temperature (K) 


FIG. 4. At left, effective barrier heights as a function of temperature in (a) diodes not annealed 
at temperatures higher than 300 °C; (b) annealed at 350 °C; (c) annealed at 450 °C; full curves 
are the best fits according respectively to Eq. ([5]) or Eq. ([9]) which results essentially in the same 
behavior for the effective barrier. The point at the highest temperature in data (b) must be 
discarded because it has been measured at the same temperature as the annealing one, involving 
an additional annealing time which has induced a further evolution of the interface towards its final 
state (c). At right, ideality factors for the same diodes after the same respective post-treatments 
in (a’), (b’) and (c’) ; the best fits to the ideality factor data obtained with the help of Eq. (fTTI) are 
displayed in full curves for (b’) and (c’) and as a dash-dot curve for (a’), but using a non physical 
parameter in the last case. The full curve adjusted on (a’) data exploits Eq. (fT7)l and related ones 
with more meaningful parameters at the expense of a worse agreement with data at the lowest 
temperature (see text). 
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realistic simulation of some potential fluctuations at interface. But within the framework of 
the simple approach developed in the present section, it is already possible to derive the two 
parameters ^b{oo, 0) and cr$(0) at zero voltage from the £t of these quantities according to 
Eq. ([5]), which are renamed respectively and for sake of simplicity and displayed in 
Table H for the three types of interface. The average barrier height at inhnite temperature 
<h^oo decreases by about 0.5 V after each anneal while the standard deviation of the 
barrier fluctuations goes through a maximum after annealing the interface at 350 °C. But 
the dehciency of the present model is conhrmed by the different behaviors of the ideality 
factor with temperature plotted in Fig. 0] which cannot be explained. These discrepancies 
conhrm the need for a more accurate model of the effect of potential fluctuations at interface 
upon current-voltage characteristics. 

Meanwhile, an unusual result which ensues from the calculation of the current densities 
necessary to draw the htting curves in Fig. [2] can be already pointed out. At some current 
threshold, flat bands are reached (see the arrows in the same hgure) and it is necessary 
to consider that accumulation occurs at the interface to explain the continuous increase of 
the current densities beyond this threshold. Such a behavior matches Eq. ([T]) which states 
the proportionality of the current to a velocity and carrier concentration at interface. For 
example in the case (c), flat bands are reached at a voltage of 0.63 V, corresponding to an 
initial band bending of 0.97 — 0.38 = 0.59V (here Vdop=0.38 V) plus the barrier increase 
of 0.04 V demonstrated in the next section. The corresponding current density calculated 
directly with the hole concentration given by formula (1) in ref.[l|, assuming a compensation 
of 20% (donor concentration Nd = lO^^cm”^) and Eq. ([T]), with a collection velocity ten 
times smaller than that for an intimate contact, yields 0.4 A/cm^, very close both to the 
experimental one and that calculated by Eq. (01). In cases (b) and (a), still at 300K, flats 
bands are reached respectively at lower current densities indicated in Fig. ([2]), which cannot 
be calculated with the help of the initial model because the apparent barrier heights deduced 
presently are largely underestimated and would produce irrelevant band bending. Therefore, 
assessing these values needs to wait for more accurate and physically relevant models which 
will be developed in the next section. In Eq. (E]) and subsequent ones, the barrier height 
has been introduced for convenience and it was often believed as an impassable limit in 
Schottky diodes. Such a statement is not true as everyone can be convinced if he realizes 
that Eq. ((2]) and subsequent ones all derive from Eq. ([1]) where the current is controlled 
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by a carrier concentration. Band curvature reversal and occurrence of accumulation regime 
are possible and even necessary to reach high currents, up to more than thousand times 
those got at flat bands, in these interfaces made on oxygen-terminated diamond, because 
hrstly of the low values of the effective Richardson constant, secondly of the wide forbidden 
bandgap of diamond and thirdly of the large ionization energy of acceptors. In such a 
regime, diffusion flow is inverted close to the interface, the barrier seen by majority carriers 
is canceled and consequently, the models relying on Eq. ([2]) and following ones, and on 
the properties of barrier height inhomogeneities, specially their voltage dependence, become 
inappropriate. Conversely, current crowding near the pinch-off patches with a lowered barrier 
may become a major source of non ideality. These mechanisms will be either incorporated 
into the forthcoming models or else, kept in mind to assess the applicability of the following 
developments. 


IV. POTENTIAL BARRIER INHOMOGENEITIES AND IDEALITY FACTOR 

A. Potential fluctuations due to randomly distributed discrete charges and dipoles 

The electrostatic potential experienced by mobile carriers originates from hxed charges 
scattered inside the space charge zone and interface. For modeling purposes, it is often 
derived from constant bulk charge concentration in the space charge zone and homoge¬ 
neous charge density at interface. This view is correct for mechanisms characterized by 
lengths larger than the Debye length Xd = \/ ^sc^o^bT /(coA^), where esc and Eq are re¬ 
spectively the relative permittivity of the semiconductor and the vacuum permittivity, and 
when average laws are evaluated, but may fail if some physical property of carriers experi¬ 
ences microscopic fluctuations. Indeed, when the effects of potential inhomogeneities are at 
stake, the discrete nature of hxed charges accommodated within the metal/semiconductor 
or metal/oxide/semiconductor (MOS) structure must be taken into account and the average 
distance between these discrete charges have to be compared to Ad. Here, Xd amounts to 
40 nm at 300 K and 52 nm at 500K. Three locations have to be considered: (i) inside the 
depletion zone, the average distance between ionized acceptors is 58 nm, very close to Ad- 
Only those which are at such comparable distances from the interface influence the potential 
barrier height, as clearly evidenced in a following subsection, (ii) The oxygen terminations 
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bonded to the carbon atoms of the (001) diamond snrface are making up the negative part 
of the interface dipole, really present at the interface up to some annealing temperature, 
an inverted situation in comparison to what happens on the hydrogenated diamond sur¬ 
face, due to the sign change of the electronegativity difference between carbon and oxygen 
or hydrogen respectively. The positive counterpart of this dipole lies within the hrst and 
second carbon layer just below the surface, resulting in a potential jump through the sur¬ 


face of Aq ddi Ndi/(eiSo) where Aq is the charge trans: 
site density and £* the relative interface permittivity 


er, ddi the dipolar distance, Ndi the 


24| . This potential jump is directly 


reflected by the electron affinity change provided by oxygen terminations. In average, it can 
be estimated close to 1.7 V because the smaller difference between electron affinities of the 
free (2x1) reconstructed surface and the oxygenated surface |33j is probably due to already 


electron depleted bonds in the last carbon layer of the former one. Taking into account the 
30 percent coverage measured by X-rays photo-electron spectroscopy [2l| on our samples 
treated by UV-ozone, equivalent to 5 x 10^^ sites per cm^, a whole distance normal to the 
interface of roughly ddi = 0.1 -|- 0.2 = 0.3 nm and Si = 2, the estimated charge transfer turns 
out to be 2x10“^° C, that is 'd'=0.13 electronic charge. However, it is known hrstly that 
there are several oxygen terminations, at least ether (C-O-C), carbonyl or ketone (C=0), 
hydroxyl (C-OH) and carboxyl (0=C-OH) groups [34143 7l| as individual adsorbates which 
can be simultaneously present on the diamond surface, and possibly epoxide (bridged oxygen 
between two adjacent carbons) 38(] in case of UV-ozone treatment as done on our samples. 
The electron affinities calculated theoretically for the three hrst terminations are different 
3a, [SOj. Secondly, all the sites are not occupied by ether or ketone terminations (in our 
case, about 30 percent 2l|, like when oxygen plasma treatment is applied to the diamond 


surface) but when adjacen 


sites are provided with an oxygen termination, polymerization 


along chains may occur 3^ . Thirdly, comparison between the experimental electron affinity 

oriented diamond surface [33j, |40|], near 1.6 eV or 1.7 eV, and 
39| (the ’’bulk” one in the sense of the Sque publication. 


of the oxygen terminated (100 
those calculated theoretically 


36 


that is the difference between the bottom of the conduction band and vacuum level at the 
surface) shows that the latter are about twice the experimental ones for ether and carbonyl 
terminations and close to zero or negative like for the hydrogen termination in the case of 
hydroxyl. It is therefore quite clear that a mixing of at least these three types of termi¬ 
nations is necessary to reconcile experimental and theoretical electron affinities as long as 
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those measured by photo-electron spectroscopy are averaged values. Therefore, the dipole 
strength resulting from these various adsorbates and able to modify the potential barrier 
height is far from being homogeneous at least over a scale length of few nanometers and 
also along a much larger scale if the various terminations are arranged in macroscopic do¬ 
mains, resulting in a non-uniform local electron affinity. An other reason which may lead to 
scattered patches with different barrier heights would be the chemical inhomogeneities like 


often present when carbide formation occurs 


13| . which might also happen if the chemical 


composition of the interfacial oxide layer or its thickness were not homogeneous throughout 
the interface. And last but not least, extrinsic defects such as emerging dislocations cannot 
be excluded as sources of locally lowered barrier height like revealed at CoSi2/n-Si(001) 


interfaces by ballistic electron energy microscopy 


41| . On the contrary, the intrinsic dipole 


due to charge transfer across the interface because of a possible misalignment of the charge 
neutrality levels in the two materials is supposed to be homogeneous because it results from 
the occupation of gap states which are delocalized in the plane of the interface and will not 
participate in potential fluctuations, (iii) As evidenced in the following, discrete charged 
centers exist inside the oxide of MIS structures and are expected to be also present in the 
thin interfacial oxide of the Schottky contacts. These cases are treated quantitatively in the 
next subsections and are brought face to face to experimental and simulation results. 


B. Evidence of charged centers in oxide layers 


To check quantitatively the presence of discrete charged centers inside the oxide, a MIS 
structure Al/Zr02/p-type (001) oriented diamond has been fabricated. The 25 nm oxide 
layer has been prepared by atomic layer deposition at 100 °C. The p-layer has a thickness 
of 2 pm, contains some 10^^ B/cm^ and has been epitaxially grown on top of a 200 nm p^^ 
diamond layer itself in epitaxial relationship on a (001) oriented Ib substrate. After etching 
the p-layer and before the deposit of the last stack Al/Zr 02 , an ohmic contact is fabricated 
around the sample on the heavily doped (p^"*") layer, ensuring a series resistance close to 
a few khl. In this way, the proper time constant of the MIS structure is inferior to 1 ps. 


Such a property allows safe measurements of the sta 
at a frequency of 1 MHz, as described in reference 


dc, pulsed and transient capacitances 


42|. Results of the reciprocal square 


capacitance as a function of voltage and the isothermal spectra of projected interface states 
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as a function of the hole emission time constant deduced from the hrst real Fourier coefficient 
of the capacitance transients at two time window durations are displayed respectively in 
Fig. El^a) and Fig. El^b) and (c). From the extrapolation of the linear part of the curve 
in Fig. El^a) to zero ordinates, a flat band voltage VpB oi —4.1 V is deduced whereas a 
value equal to <Fms = — Xsc ~ EQ/2eQ — {Ei — Ep) /cq = —2.4 V would be calculated 


with the help of the electron affinity model and = 4.28 V for the A1 work func 


43|, Xsc = +1.6eV for the electron affinity of oxygen-terminated p-type diamond 


33 


ion 


m 


Eg = 5.45 eV for the diamond band gap and 0^ = {Ei — Ep) /cq = 2.35 V deduced from the 
difference between intrinsic and Fermi energy levels in bulk boron-doped diamond, assuming 
no charge in the oxide. The discrepancy comes from this last hypothesis which is not verihed, 
implying that the real situation comprises a positive charge inside the parenthesis of the last 
term in VpB = ^ms— {Qis{0) + Qox) /Cox-, where Cox is the oxide capacitance density, Qis(O) 


and Qox being respectively the charge per unit area in interface states at flat bands and the 
product of the total charge per unit area within the insulator by the ratio of the barycenter 
of this charge to the insulator thickness 42(]. To reconcile the two values, a total number of 
positive charges per unit area must amount to 2x10^^ cm“^ within the 25 nm thick oxide 
layer. 

The centers which accommodate charge may experience change of their charge state 
if some mechanisms are occurring, like carrier emission from interface states lying in the 
common forbidden gap of the two materials, and also transport by hopping and tunneling 
across some thickness of the insulator. Such effects give rise to transient capacitances in the 
MIS structure because of the rearrangement of all the individual electric dipoles necessary 


to comply with t 


le rule of a constant whole electric dipole if the MIS is submitted to a 


constant voltage [4^. Isothermal spectra are displayed in Fig. |5](b) and (c) as a function 


of the characteristic emission time. They are obtained from the Fourier transform of the 
transients which are recorded after two Filing pulses lasting respectively 1 ms and 100 ms. 
They show a two peaks structure, the slower peak being higher for the longer Filing pulse. 
This response can be viewed as the projected spectrum of states on the interface energy scale 
inasmuch both states located at interface and some oxide states as well may exchange holes 
with the valence band of diamond. DiFerential spectroscopy is here unnecessary just because 
emission from discrete interface states in the common bandgap largely dominates that from 
a hypothetical continuum since peaked structures clearly appear in Fig. [5l and an accurate 
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FIG. 5. (a) Square of the inverse capacitance density as a function of the gate voltage ; (b) and (c) 
curves are isothermal transient capacitance spectra respectively after a 1 ms and a 100 ms filling 
pulse, recorded at a gate voltage of 3 V and temperature of 362 K. 

calibration of the energy scale is not needed. The smaller hlling pulse duration tp = 1 ms 
has been chosen to ensure complete hlling of interface states even with capture cross sections 
as low as 10“^® cm^, whereas the larger tp = 100 ms is able to involve additionally either 
injection or extraction of carriers trapped in discrete states localized within the insulator at 
some distance from the interface. The important point consists in the different behaviors 
followed by the two peaks as a function of tp-. the hrst occurring at the shorter emission 
time is very weakly sensitive to tp and consequently related to an interface state whereas 
the second one increases with tp, involving carrier transport through some insulator slice. 
This is the demonstration that at least one discrete state exists inside the oxide. An order 
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of magnitude of the distance d of these centers from the interface can be estimated from 
the inverse of the attempt-to-escape time checked by Jain and Dahlke for traps in Si/Si02 
interfaces H and given by z/q exp(—gd), with z/q ~ 10^® s ^ and a ~ l.OSA yielding 
roughly d ~ 3 nm. The reciprocal attempt-to-escape time at zero insulator thickness uq 
and the quantum transmission coefficient q are conditioned respectively by the quantum 
uncertainty about the potential well or barrier height, also close to 2h/G where G is the 
line width of the resonance and tunneling transmission coefficient which depends on 
the square root of the previous quantities. Because the relevant energies are still about 
1 eV, the order of magnitude of uq and q are expected not to be largely different in the 
present case. The evaluated distance d is close to one tenth of the whole insulating layer 
thickness, a ratio similar to that of the center densities deduced by the transient and static 
capacitance methods in the MIS structure. Consequently, these point defects behaving as 
repulsive centers for holes and located inside the oxide layer at the Schottky interfaces under 
study might not be ignored unless other charged centers or dipoles prevail. 


C. Effects of local barrier enhancement due to random elementary charges and 
dipoles close to the interface 


Electrostatic potential along x and y directions, parallel to the interface, and z, normal to 
the interface, can be calculated by adding three terms as a function of normalized distances 
^ = x/Xd, y = vIXd and C, = z/Xd^ and band curvature U in volt: the average potential 
of the depletion zone (y\/U — ksT /cq — (^/ksT /(2eo) j ; the image force potential —UimIC 
with Uim = eo/{16nesc^o^D) and the potential due to either repulsive elementary charges 
positioned symmetrically at .^ = ±.^ 0 ) y = 0, and beyond the interface plane by a distance 
Co, that is 




+ 


4m,; 


y(5 - + (c+ Co)'' 

or replacing either one or two of the previous terms by 


+ (C + Co)^ + (C - Co)2 

at either one location -|-Co or both locations ±Co to simulate either one or two elementary 
dipole due to negative oxygen termination beyond the interface and its positive counterpart 
on the diamond side. A 3-dimensional view of potential barriers in several cases is given in 
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FIG. 6. Sketch of the potential barrier experienced by majority carriers with (a) two dipoles across 
the interface with a repulsive elementary charge on the semiconductor side; (b) one dipole of the 
same kind across the interface and one repulsive elementary charge beyond it; (c) two repulsive 
elementary charges beyond the interface; (d) randomly scattered attractive elementary charges 
inside the depletion zone, inducing a saddle point of the potential energy, visible as a lowered 
profile of the top of the barrier near the potential well which is the closest to the summit. 


Fig. El the last one being obtained with an inverted sign before (o and four random locations 
of elementary charges within the depletion zone. The three hrst cases are respectively 
devised to mimic the mechanisms assumed to be at stake in the transport of carriers at 
interface, resulting from the repulsive effect of charges located just at interface or close to 
it inside diamond: (a) with negative oxygen terminations inducing interface dipoles which 


1 / 
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are supposed to dominate; (b) with a mixed influence of repulsive oxide charges and dipoles 
similar to those of the previous case; (c) with only repulsive oxide charges. In these three 
cases, elementary charges or dipoles enhance the barrier locally. On the contrary, in the 
fourth case (d), barrier lowering 5$™'^ may occur with a saddle point at normalized distance 
Csad from the interface which can be calculated exactly from the solutions of a quartic 
equation resulting from the cancellation of the potential derivative along the ( direction. 
However, a rough estimate, easier to derive, relies on the geometric mean of the two extreme 
solutions of the quartic equation, and can be expressed as (sad ~ \/x{U) [Co — 2x([/)] with 
xiU) = ‘ results in j 4 (Co - Cad)- The 

resulting effects of this barrier height lowering will not be treated in detail here. In the 
three first cases, the distribution function of potential barrier heights is then calculated with 
a fixed value of Co and random values of Co around a mean Co according to a normal law 
with a standard deviation cr^, to assess the possible values of ^B(T,Vj) and cr$(V^), the 
characteristic parameters of the inhomogeneous barrier. 


In a first step, the parameter cr$(0) is adjusted to match the values displayed in Table |T] 
at a medium temperature, typically 500 K with realistic values of (o, Co and a^. More 
specifically, Co-^d spans the range 0.2 to 0.46 nm and (Co-^d)”^ the range 1.5 x 10^^ cm“^ 
to 3 X 10^^ cm“^, the two later limits featuring respectively a positive charge density inside 
the oxide layer and the dipole density at the oxygen terminated surface. Accordingly, the 
average distances Co between point charges or dipoles span a range limited by 0.08 and 0.007 
(in Xb units) respectively. Therefore, the barrier height fluctuations are occurring along 
a microscopic scale, well below the Debye length. In a second step, the band bending is 
changed and the values of ^^(T, 77) and cr$(f/) are evaluated as a function of the band 
bending U. In all the cases and whatever the temperature, as suggested in |29[, the increase 
of the average barrier height ^B{T,Vj) and the decrease of the standard deviation a^iVj) 
with increasing forward voltage and correlative decreasing band bending is fully confirmed, 
although the barrier height distribution is not symmetric and cannot be strictly depicted by 
a normal law. The trends followed by the two terms of Eq. ([5]) both contribute to an increase 
of the barrier height with forward voltage which can be also interpreted as an increase of the 
ideality factor, thus inducing a lower current density in comparison to that anticipated with 
the help of a constant ideality factor like sketched with dash-dot lines in Fig. [2l At variance 


with the hypothesis used in ref. 


29| . neither the variation of <hs(T, Vj) nor that of aq.iVj)" 
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FIG. 7. Ideality factors deduced from the microscopic barrier inhomogeneities model for the three 
types of interface at 300 K (triangles), 400 K (squares) and 500 K (circles) as a function of the 
band bending U. Curve (i) is the best fit of data points by 1 + 1/P{U) with P{U) = 136.81/ — 
76.41/^ + 28.417^. The curve (ii) displays the barrier height variation deduced from (i) with help 
of Eq. ([6|) which implies that / dU = —1/ [1 + P{U)]. Dash-dot curve (iii) shows the image 
force lowering variations. 

is linear with Vj or U, a law which would have resulted in a constant ideality factor, in 
contradiction both with experimental and simulation results. An other major feature is the 
relative insensitivity of the ideality factor to the density of point charges or dipoles (^o^d)”^, 
the former starting from 1.01-1.02 for band bending corresponding to zero voltage bias and 
growing up to 1.05-1.15 when the band bending U becomes close to a few ksT/eo like shown 
in Fig. [TJ irrespective of the temperature. The barrier height variations which are deduced 
from the ideality factor match very well those of the image force lowering as clear in Fig. [7] 
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too. This result coming from the simulated inhomogeneous barrier height with repulsive 
elementary charges or dipoles demonstrates that the image force lowering and its changes 
with the band bending is the main source of the variation of the average barrier height with 
the applied voltage, because barrier height enhancement does not contribute signihcantly to 
the whole barrier experienced by carriers responsible for the current, except for a reduced 
effective area. Such a result validates this hypothesis already done in ref. j^. The whole 


variation amounts to about 40 mV like inferred in ref. 


24| for the image force lowering alone. 


It is thus possible to derive a universal law featured by a solid line in Fig. [7] that may be 
used to derive current density-voltage characteristics with the help of Eq. (|1]) and ([H]). The 
hrst important result coming out from this simulation turns out to be the better agreement 
of the simulated current density sketched in solid line in the inset of Fig. |2]in the case (c) 
where the barrier height correction has been taken into account. The second one is the 
insensitivity of the ideality factor to temperature, a property which matches rather well the 
behavior of the interface (c) in Fig. |4](c’), in contrast to the absolute variations in data (a’) 
and (b’), although the ideality factor (c’) remains slightly higher than the values evaluated 
from the image force lowering alone. Taking all the arguments evidenced and discussed in 
the present subsection, one must conclude that microscopic barrier inhomogeneities resulting 
from repulsive elementary point charges or dipoles which have been however, experimentally 
evidenced, cannot account for most of the properties of the interfaces presently studied, and 
that an other type of barrier inhomogeneity must be present. 


D. Barrier height and ideality factor of interfaces with lowered barrier patches of 
size either comparable to or smaller than the Debye length 


The calculation of barrier height fluctuations due to barrier height lowering Ap at circular 
patches with radius Rq or linear ones with width Lq and their effect on the current-voltage 
characteristics relies on the same principle as previously, except for the weight of each dipole 
which has to be 27iesc^oApRQ to induce the desired barrier change at circular patches which 
will be only considered in the following because close to those happening naturally in inho¬ 
mogeneous interfaces. Additionally, only patches with a lowered barrier are considered as 
demonstrated in the previous subsection. This method, known as the dipole-layer approach, 
has been implemented analytically by R. T. Tung 30| and summarized later by W. Monch 
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24j |. For a single circular patch with radius Rq and barrier lowering Ap embedded into 
a homogeneous area characterized by a constant barrier height a single parameter 

containing the two characteristic quantities Rq and Ap, which is called the patch param¬ 
eter 7 p = 3(Ap can be dehned. Integrating Eq.([T]) over the whole area A with 

the hypothesis of a parabolic shape of the potential around the saddle point, which is the 
lowest potential of the patch, restores the form of Eq. ([3]). But the new barrier height is 
jjiinus an effective barrier lowering at the potential saddle point, the ideality 

factor n is replaced by 1 and an effective area A^ff = dvrAf^F/d affects the current expres¬ 
sion. The weight of the patch is proportional to the dimensionless parameter F = ^ 31 J 2/3 

where rj = esc^o / (eo-A^yi)- Both the effective barrier lowering 3T U = jplU / 7 ]^^^ and area 
Aeff = [v / are dependent on the band bending U = — Vj — Vdop, while 

the ideality factor is approximately 1 + F, therefore increasing with the forward voltage. 
The hnal step consists in introducing a distribution of a variable patch parameter. In the 
Tung’s original work |30|, two cases were worked out : (i) a sharp distribution of lowered 
barrier patches approximated by a single patch parameter 70 and (ii) a broad distribution 
featured by half a Gaussian law which limited the 7 ^ range to the values pertaining to 
barrier height lowering alone, that are 7 ^ ^ 0. This choice was justihed by the negligible 
contribution of patches characterized by a barrier height enhancement to the whole current 
as we checked in the previous subsection. In order to generalize such an approach, allowing 
to match continuously any intermediate behavior between the two previous cases, the patch 
parameter distribution function will be chosen as N ( 7 ^) = character¬ 

ized by a standard deviation ap of the patch parameter around the central value 70 and an 
areal density Dp. Depending on 70 , the function N^'jp) depicts either a peaked distribution 
around the central patch parameter, or half a Gaussian law decreasing from its maximum 
at 7 p = 0 when 70 = 0 or even a decreasing Gaussian tail with a non zero initial slope if 
7 o < 0. After summation of all the currents dependent on any positive value of 7 ^ from 
zero to inhnite, this procedure further allows to derive the following expression where only 
the last factor inside braces contains the effects of inhomogeneities while the others give the 
ideal thermionic current density of a homogeneous barrier height 
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Bnt the cnrrent density deviates notably from the ideal thermionic law only if the factor 
inside the braces of Eq. ([7]) is larger than nnity. As one is interested in this case, a reasonable 
approximation consists in neglecting the nnity terms inside the braces, so that the effective 
barrier height can be dednced from the exponent of the exponential factors, which turns 


out to be i 


ep ^ X^sad ^2 _ 7 o_ 


inside the braces of Eq. ([7]) and logarithm of the pre¬ 


exponential factor : 
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This last equation is similar to Eq. (|5]) with the new quadratic average <5$^ > of barrier 
fluctuations 
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and an additional term which has a little influence due to the logarithm and a null value 
for Dp in the 10 ^° cm“^ range, but may give a non negligible contribution to the ideality 
factor appearing in next formula when the barrier inhomogeneity level is low. But the main 
issue is the dependence of <6 ^b> on temperature which is different from that given by 
Eq. ([5]) and results in a shift of the high temperature asymptote of the effective barrier in 
comparison to by a quantity proportional to 70 , leading to an apparent extrapolated 
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From Eq. 


and ([9]), the ideality factor can be calculated as: 
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where is between 1/3 and 2/3 according to the respective dominance of either the hrst or 
the second term in < <5$^ >. This ideality factor increases with decreasing band bending 
and therefore with increasing applied voltage as long as flat bands are not reached, so that 
current-voltage characteristics plotted in a semilogarithmic diagram experience a downwards 
curvature, as in Fig. [2] (a) and (b). This model ceases to be valid neat flat bands and beyond 
because the current calculation relies only on the carrier flow at the saddle points of the 
potential barrier, neglecting the areas where the barrier is or higher. It must be also 
stressed that these ideality factor and effective barrier height must be evaluated around the 
same band bending U = Ub if comparison at different temperatures or in different diodes 
are made. This is often done by assessing the corresponding ideality factor Ub at the same 
current density although such a rule is almost true only at vanishing band bending. It is 
important to notice that various behaviors of the ideality factor as a function of temperature 
can result from Eq. mil. with a trend toward a constant value when the second term inside 
the brackets dominates at elevated temperatures as one can see in Fig. ([SJd’), the last term 
in Eq. (1TT]1 being negligible except near flat bands and for ideality factors below 1.1. Either 
neglecting the last term in Eq. (I9]) or lumping it into a relationship between the 

effective barrier height and ideality factor nb at the same band bending Ub for all the diodes 
investigated at a given temperature among a set of diodes, comes from Eq. (ED, (IinD and 
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Therefore, the difference between the effective barrier height and a quantity shifted 
from the homogeneous barrier mainly by a term proportional to the central patch 

parameter tq, follows a linear relationship regarding the ideality function (n^ — l)/nft. It 
does not depend on the standard deviation cXp of the patch distribution, except for the last 
term including the factor ki^, which is roughly constant for a given patch distribution. At 
constant temperature, this last term is indistinguishable from and adds a minor error 
to it. Conversely, the term proportional to To may contribute to an important shift of 
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Sample type and 

<^-PP 

(Tp 

70 

^hom 

Dp 


annealing temp. 

V 

Vl/3 

Vi/3 

V 

_2 

cm 

V 

(a) < 300 °C 

2.04±0.14 

(9±3)x10-5 

(1.6±0.4)xl0-^ 

2.26±0.16 

(3±l)x 10® 

2.38±0.02 

(b) 350 °C 

1.68±0.1 

(9.0±1)x10-5 

(1.85±0.3)xl0-^ 

1.87±0.1 

(6±2)x109 

1.90±0.02 

(c) 450 °C 

0.93±0.04 

(5.5±1.5)x10-5 

(1.3±1.2)xl0"® 

0.94±0.04 

< 10^ 

0.99±0.01 


TABLE II. Parameters of the Schottky barrier in the three interfaces according to the models 
detailed in section lYl with and ap deduced from the adjustment of barrier height as a 

function of temperature according to Eq. Q. Then, 70, and are derived from ideality factors 
fitted as a function of temperature with respective help of Eq. m in the case of interfaces (b) and 
(c), and Eq. (I16bp for interface of kind (a), while Dp and come from the direct adjustment 

of the current densities shown in Fig. ([2|) respectively to Eq. ([7|) for interfaces of kind (b) and (c) 
and Eq. (jl4ll for interface of kind (a) as a function of voltage. 


the apparent extrapolated barrier at unity ideality factor in comparison to since this 

/ \ 1/3 

Y f — j and this correction might not 


\hom 


apparent barrier turns out to be 
be ignored in general. It can also introduce some dispersion of the data when a family of 
diodes is studied. The effective barrier height is approximately linear with rih only for 
typical values smaller than 1.4 like noticed in ref. 2^, so that it is better to make plots of 
the apparent barriers as a function of {rib — for sake of generality, as done in Fig. (jH]), 
which will be commented in the next subsection. However, it must be stressed that all the 
expressions derived for the effective barrier height and ideality factor within the framework 
of the present model, that are Eq.(l9]) to ([12]), assumed that the factor inside the braces of 
Eq. (|7]) is largely greater than one. Such a condition implies that the patch density Dp 
must be typically much larger than 10^ cm“^, so that one has to wonder whether the value 
deduced from adjustments has really a physical meaning in each situation. 

In one diode alone, the hngerprint of the kind of barrier inhomogeneities breaks out in 
the behaviors of the ideality factor and to a smaller extent of the effective barrier height as 
a function of temperature, like shown in Fig. (jlj), where evolutions of the ideality factor are 
markedly different for the three interfaces in Fig. dlti’)) dlb’) and (jH::’). All these features can 
be reproduced by the laws demonstrated above, specially the ideality factor variations with 
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FIG. 8. Effective barrier height as a function of {uf, — l)/nh, where is the ideality factor at band 
bending 17;, for 95 diodes fabricated on three distinct diamond stacks elaborated like described in 
section [ni respectively identified by squares (I), pentagons (II) and circles (III). Samples I and II 
have not been intentionaly annealed while sample III has been annealed at 350 °C. Dash lines are 
best linear fits of grouped data lA, IB, IIA, IIB and III according to Eq. (fT^ . 


temperature which are strongly influenced by the central patch parameter 70 in Eq. (ITTH . 
In the case (a), a good agreement with Eq. ffTTj) can be obtained as depicted by the dash- 
dotted curve in Fig. (jl^’) at the expense of a negative central patch parameter which would 
induce an apparent barrier larger than the homogeneous one obviously leading 

to an incoherent physical situation. Moreover, the direct adjustment of the current density 
according to Eq. ([7]) turns out to be impossible because this expression is unable to account 
for the curvature of the current density-voltage characteristic (a) in Fig. ([2]) below flat 
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band voltage. An alternative model able to involve a more appropriate physical picture is 
discussed below for interfaces of kind (a). Data collected in 95 diodes prepared on three 
distinct diamond stacks which underwent the same pre-treatment before metal deposition 
and displayed in Fig. (E]), give a hrst indication: if so different groups of data are possible for 
diodes whose diameter is 100 pm, specially the IB against lA and IIA, related to interfaces of 
kind (a), barrier heights extrapolated at unity ideality factor and averaged on such area must 
be different. Consequently, the lateral scale of barrier inhomogeneities may no longer be the 
Debye length but much more, and the approximations of a parabolic potential at the bottom 
of the lowered barrier patches depicted by the patch parameter 7^ would cease to be valid. 
Secondly, strong variations of the ideality factor with temperature like in Fig. suggests 
that diamond resistivity may play a role because acceptors are progressively ionizing in this 
temperature range. 


E. Barrier domains with multiple sizes and inhomogeneities of the current density 

When the lateral size of barrier inhomogeneities largely exceeds the Debye length, the 
previous model relying on a parabolic potential at the bottom of the patches with a lowered 
barrier is no longer relevant. These different areas with distinct barrier height can be treated 
as independent of the band curvature because the pinch-off of the channels where the current 
flows becomes negligible. Then, the total current is simply obtained by summation of all the 
currents like in a parallel circuit of the corresponding paths. But such a model is unable to 
£t the downward curvature experienced by the logarithm of the current density as a function 
of voltage in interfaces of kind (a) and yields merely an ideal behavior as soon as the lowest 
barrier is shifted below the others by some ksT/eo because the corresponding current flow 
dominates the others. In fact, this current inhomogeneity due to carrier crowding at the 
saddle points and lowest barrier domains remains when the ^ coordinate move away from the 
interface and therefore have a dramatic influence on the total voltage applied to the junction 
for a given current, specially when carrier freeze-in occurs, inducing large variations of the 
ideality factor with temperature. This effect is relevant both for barrier domains much larger 
than the Debye length where no pinch-off occurs and for pinched-off patches, and the voltage 
drops in these shrunk channels must be subtracted from the total applied voltage. In the 
pinched-off case, the diameter of patches is always much smaller than the lightly doped layer 
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thickness ti so that the spreading resistance can be approximated by — j £=— with a good 


^/A 


If,’’-' 


accuracy 4^, p{T) being the resistivity of the lightly doped diamond layer. Therefore, the 
non pinched-off and pinched-off forward currents which flow respectively through individual 
areas Ai and A^Jj ^ under an applied voltage 14 can be derived from the generic thermoionic 
law but with band bending equal to the homogenous potential barrier minus the relevant 
barrier lowering and applied voltage 14 from which the ohmic voltage drop along either non 
pinched-off or pinched-off channels inside the depletion zone is subtracted: 
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the lowest barrier spreading over the areas Ai with size much larger than the Debye 
length being — A<f). If the current inhomogeneity remains the same along the whole 

depth of the lightly doped layer like assumed in ref. Ml], the voltage drop along the 4^ non 
pinched-off channel is simply proportional to a resistance p(T)/Aj where Aj = Ai/ti. This 
hypothesis is true in the limit of a diameter of the area with a lowered barrier height well 
larger than ti, a condition equivalent to \/Aijn S> fz/2. But because this diameter can span 
a large range, from values well below ti jm to much wider ones, the corresponding spreading 
resistance should take the general form 46] arctan(2f// A^j'n') and A, should become 

^yJ-KA^j arctan(2b/ yjAij'^'). It must be noticed that the last term in the exponent of the 
expression of If and in f/j is added to 4)^°"* and therefore plays the same role as a negative 
central patch parameter jq in the effective barrier experienced by the whole current flowing 
through the pinched-off patches and given in Eq. ([9]). This fact explains why neglecting the 
conhnement of current in channels which stretch out along some part of the lightly doped 
layer thickness had to be compensated by the occurrence of a negative hgure for 70 in the 
simulation results of the previous subsection, leading to an incorrect value because of an 
inappropriate model applied to interfaces of kind (a). Such a fictitious negative figure for 70 
was necessary to obtain an acceptable, but misleading, agreement of the experimental data 
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with the model developed in the previous subsection since firstly, the resulting value 
larger than would contradict the hypothesis which assumed lowered barrier patches 

and secondly, the same model was unable to match the curvature of the current-voltage 
characteristics. 

Averaging all the currents over the whole diode area Ad is not tractable analytically. A 
reasonable approximation consists in calculating separately the total current and an estimate 
of the whole conductance within the shrunk channels from the integration of all the parallel 


conductances Aj/p(T) or AF^^j J-k/ p{T) appearing in the last term of the exponent of 


each expression (1 13 ah or fll3bl) . weighted by the distribution Nip^p) = 


Dr, 


p‘lixa 


■ exp 


(7p-7o)" 

2^2 


in the later case. It is also assumed like previously that only the lowest barrier for the current 
flow is efficient and that the barrier lowering is approximated by its quadratic average. For 
non pinched-off channels, the total conductance would result from the integral of Aj/p(T) 
weighted by an unknown distribution, so that it is more convenient to introduce a hctitious 
area allowing to write the corresponding voltage drop ^y^-on-p-o^ 

of the pinched-off currents, two integrals must be evaluated numerically as a function of 
7o/(\/2ap). Their ratio appears as a fraction in the next formula and gives the exact value 
within an error smaller than 10 % . With = 'Yhi representing the true area 

of the non pinched-off zones, the average current densities in the whole diode area become 
respectively for the non pinched-off and pinched-off channels: 
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It must be noticed that the patch density Dp does not appear in AVp~° because both 
the conductance and current are proportional to it, so that it is canceled in the frame¬ 
work of the previous approximation. From these expressions, it is easy to guess that the 
lower the temperature, the higher the effect of the smallest barrier, like in previous mod¬ 
els. If the non pinched-off current dominates, the apparent barrier lowering would be 
A$ — ^^\n(Ad /inducing a linear decrease of the effective barrier as temper¬ 
ature is lowered, whereas alternatively Eq. (|9]) would be valid, at least when AV remains 
small regarding I 4 . At sufficiently high temperature, the increase of the effective barrier 
would be linear with temperature when the non pinch-off current dominates, an effect which 
is not clearly detected in Fig. indicating that the pinch-off current might stay the 

major one. Such a view is conhrmed by the absence of any ”S-shaped” current-voltage 
characteristic in Fig. (E^), which otherwise would reveal the existence of two parallel non 
pinched-off channels like happening in some Ti/4Il-SiC diodes j^, each of ones having its 
own barrier height and series resistance with the corresponding condition A^ji: S> tz/2 
indeed achieved. 

From experimental data displayed in Fig. ([H]), it seems that the general trend involving a 
linear relationship between effective barrier height and ideality function is still valid. Such 
a property can be checked from the integration of the central expression in Eq. ([6]) where 
the total applied voltage 14 replaces the junction voltage Vj since the current density is now 
dependent on I 4 if one relies on the last expressions fll4p . It results in a linear relationship 
between the effective barrier and ideality factor function if the ideality factor does not depend 
on the band bending Ub- Such a condition is only approximately achieved, but better and 
better when higher and higher band bending are used, like those deduced from these data 
and reported in Table HTTl because the voltage drop AV remains negligible in this range. It 
is therefore possible to hnd a range of voltage and corresponding band bending extending 
over an interval Ub/K {K being a unitless constant) where the ideality factor is constant 
around Ub and integration of the second and fourth members of Eq. (El) gives the apparent 
barrier height in non ideal interfaces: 

tTiT, Ut) = - h ^ (15) 

where n(T, Ub) is the ideality factor at temperature T, around the peculiar band bend- 
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ing Ub- is the barrier height at unity ideality factor, which matches accurately 

/ \ 1/3 

III \ fQj, pure pinched-off barrier inhomogeneities as demonstrated in sub¬ 
section IIV PI in formula flT^ where the last term can be neglected, and is only close to 
without a precise knowledge of the deviation in other cases. Therefore, the observed linear 
correlations in Fig. [Hlcan still be justihed. When the voltage drop inside the lightly doped 
layer is not negligible, approximate expressions of the ideality function Y can be drawn sep¬ 
arately from the previous equations when either or Ip~° = AdJz~° 

dominates, which turns out to be respectively with either the ideality factor UaiT, I 4 ) under 
applied voltage 14 or the ideality factor nb(T, I 4 ): 
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When the ideality factor has to be studied as a function of temperature, measurements 
should be done at constant current density, because in such a case the product of aA**T^ 
and exponential factors in Eq. (1 16 ah and in the term of Eq. fll6bl) can be considered 

as almost invariant, making any parameter adjustment much easier. In each case, these ex¬ 
pressions shows that the ideality factor may become very sensitive to temperature, specially 
because it follows the variations of the resistivity p{T) which are very important in the range 
where freezing out of carrier occurs, typically from 300 K to 550 K in this study. Involvement 
of the voltage drop due to the flow of current through the narrow bottlenecks which lengthen 
the lowered barrier regions inwards the lightly doped layer is able to explain qualitatively 
the strong negative slope of the ideality factor in Fig. (jlK’)- The general case relies on the 
combination of the two kinds of current with a proportion P = P~°/+ jnon-p-o'^ 
on Eq. ([6]), giving the global ideality function: 


Y = 


n{T,Va) 


= P YP-° + (1 - P) Y 


non—p—o 


(17) 


The best match to the four data at highest temperatures is obtained with P = 72% and 
plotted as a full curve in Fig. (jl^i’)) conhrming the previous hypothesis about the dominance 
of the pinched-off current. The discrepancy at the lowest temperatures is commented in the 


37 


















following. But it must be stressed that the hrst term in Eq. (IT7|) is unable to account for 
the ideality factor variations, which are mainly due to the second term despite its smaller 
weight. It is thus conhrmed that the non pinched-off currents contribution is necessary 
to explain both the current density-voltage curvature and the ideality function shape with 
temperature in interfaces of kind (a). The more important voltage drop in these 

interfaces is also able to explain why flat bands are reached at a lower current density 
(see Fig. |2]) in comparison to interfaces of kind (b) and (c) in which the patch and channel 
densities are so high (see Table HTj) that the voltage drop remains negligible in this bias range. 


At sufficiently high forward bias, band curvature can be inverted in some areas or even in 
the whole diode, as the reader can realized from the comparison of (in Table UTll minus 
Vdop = 0.38 V to the applied voltage in Fig. (|2]), so that the previous expressions can no longer 
be used. In other words, energy bands bend upward near the interface when the flat band 
voltage is exceeded. From Eq. (HB, barrier lowering vanishes at flat bands for P ° while it 
stays constant for and one can think that the effective barrier height remains more 

or less at the value obtained at flat bands beyond this threshold. But previous models which 
rely only on barrier lowering are unable to take into account the progressive disappearance 
of current inhomogeneities in the lightly doped layer due to the general increase of the 
carrier concentration at the whole interface. Such a behavior would hnally lead to a more 
homogeneous current density and a voltage drop within the lightly doped layer expected 
to be simply the product of a series resistance and total current. But even this view fails 
because voltage drop turns out to be lower in Fig. ([2]) as if the resistivity were smaller 
than the predicted one. This discrepancy is neither an artifact of the adjustment procedure 
nor a self-heating effect which would hardly be justihed because of the low dissipated power 
(typically one to some tens of mW in the range of interest), but is rather due to the injection 
of carrier by the heavily doped layer, involving a conductivity modulation due to the increase 
of carrier concentration above the equilibrium one and occurrence of a diffusion current 
starting from this heavily doped layer. This effect is evidenced from the comparison of the 
current-density data and curves simulated with a constant series resistance in Fig. ([2]) and 
it can also be checked because of the lower value of the current-density voltage derivative 
compared to the theoretical specihc resistance (not shown) at sufficiently high current. A 


detailed study relying both on experiment and simulation is published elsewhere 


48] and 
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Interface type, sample 

(V) 

Ub (V) 


and diodes set identifiers 



corrected for 70 

(a) I A 

2.27 ± 0.12 

1.05 ± 0.28 


(a) I B 

2.54 ± 0.14 

1.30 ± 0.24 


(a) II A 

2.21 ± 0.04 

0.97 ± 0.04 


(a) II B 

1.61 ± 0.06 



(b) III 

1.94 ± 0.04 

0.73 ± 0.08 

2.04 ± 0.04 


TABLE III. Schottky barriers extrapolated at unity ideality factor at band bending Ub in different 
sets of diodes fabricated on three distinct diamond stacks numbered I, II and III, according to best 
fits of Eq. (|12p and after the 70 correction for the interface of kind (b). 


demonstrates that the injection effect is more and more important as the temperature is 
lowered below the freeze out of carriers and the active layer thickness is smaller than 2 fim, 
two conditions which are fulhlled here. It had been probably observed by other authors in 
Ni/SiC diodes j^, although they did not give such an interpretation. This carrier injection 
which decreases the resistivity is still efficient at lower current density, so that the hole 
concentration at interface is underestimated even at the bias voltage where ideality factors 
are measured, at least when complete ionization of acceptors is not achieved. Correlatively, 
the effective resistivity is overestimated when its equilibrium value is taken into account 
as done in the present calculation, increasing thoughtlessly the simulated voltage drop and 
ideality factor in Fig. at the two lowest temperatures. This injection effect probably 
occurs because the lightly doped layer is only 1.4 /rm thick and would disappear for the 
larger thicknesses suited to higher voltage operation. In any cases, it must be noticed that 
all the previous models discussed and built in this work are applicable to a range of bias 
voltages and currents limited by the occurrence of flat bands in the junction but that high 
injection effects may be present near flat bands and beyond. 
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F. Characteristics and origin of barrier height inhomogeneities in the three types 


of interface 


The linear relationship between the effective barrier height as a function of the ideality 
factor function {rih — l)/nh is actually checked for diodes prepared with the same procedure 
on three different substrates either without intentional post-annealing or annealed near 
350°C in Fig. ([H]). The correction of the constant term in Eq. flT^ is done only in the case of 
indisputable relevance of the model of subsection IIVDI and as it is not certain for interfaces 
of kind (a), extrapolated barrier heights and will be considered as the same in 
table m It is noteworthy that the homogeneous barrier height can be different in two 
families of diodes (lA and IB; IIA and IIB) on the same diamond stack while it turns out 
to be very similar for two families of diodes (lA and IIA) belonging to two distinct diamond 
samples, whose the initial surfaces are all oxygen-terminated. Such results conhrm that 
oxygen terminations, which are responsible for the enhanced electron affinity of the oxygen- 


terminated diamond surfaces and highest barrier at interfaces build on sue 
expected to be organized in domains displaying distinct electron affinities 


1 surfaces, are 
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resulting in distinct barrier heights, and that the same kind of oxygen terminations can 
exist on different samples. This dispersion matches that of the various molecules and bonds 
presented in subsection IIV Al and probably also the lack of full coverage by oxygen related 
molecules. The linear relationship of Eq. flT^ or flT5|l seems to apply irrespective of the 
scale of the barrier inhomogeneities in comparison to the Debye length provided ideality 
factor is measured at high enough band bending, typically two third of the barrier height, as 
justified in the previous subsection. But because the interface annealed at 350°C (Fig. (IHlIII)) 
is a matter of the lowered barrier pinched-off patch model detailed in subsection IIVD( the 
correction which gives the homogeneous barrier from the apparent one is only exerted in this 
case in table HTTl where the parameters are collected and displayed for all the diodes families. 
It must be noticed that the coherence of the various values appearing in Tables [TTl and 
m is reasonable and that an average value of 2.04±0.04 V can be inferred for the interface 
of kind (b), only about 0.2-0.5 V below the values at stake for the interfaces of kind (a). The 
relevance of the model relying for interfaces of kind (b) on lowered barrier patches with sizes 
comparable or smaller than the Debye length and depicted by a patch parameter distribution 
centered on a peculiar value 70 > 0 , with a much higher patch density than for unannealed 
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interfaces, demonstrates that the main barrier inhomogeneities take place now along this 
microscopic scale, in contrast to what happened in unannealed interfaces. One must remind 
that the patch parameters 7 ^, ap and 70 depends on the one third power of the product 
of the square radius Rq and initial barrier lowering Ap so that similar values of the two 
last quantities will not involve necessarily the same Rq and Ap, which can move in opposite 
directions. From results given in Table [TTl such statements suggest that the 350°C anneal 
resulted in the partial and random cancellation of the dipoles 0 “^'— 0 +*^' whose stability 
is mainly governed by the enthalpy of the thermodynamical chemistry of these species and 
temperature, leaving the interface in an intermediate state and inducing new barrier lowering 
with a magnitude equal to the difference of the homogeneous potential barrier in interfaces 
of kind (b) and (c), i.e. Ap 1.1 V at a microscopic scale, and increasing considerably 
the patch density. On the contrary, the initial inhomogeneities in unannealed interfaces of 
kind (a) likely resulted from the oxygen terminations dispersion over areas with larger size 
but smaller barrier lowering and much lighter density, with typical Ap ^ 0.3 V, similar to 
the difference in potential barriers reported for various interfaces of type (a) in Table IIIIl 
To comply with the similarity of cXp and 70 for interfaces of kinds (a) and (b) in Table [TTl 
patches area has to be decreased by a factor of nearly four when going from case (a) to 
case (b), leading to a very good relevance of the model developed in subsection IIV D[ In the 
case of interfaces of kind (c) annealed at 450 °C, these dipoles have completely disappeared, 
involving a barrier height as low as 0.96±0.04 V (average from Table ITTl) and very low barrier 
fluctuations, as indicated by a central patch parameter 70 one decade smaller than that of 
other cases, which induced an almost ideal behavior of the current-voltage characteristics. 
In this last case, residual barrier inhomogeneities could be due to the random influence 
of the ionized acceptors close enough to the interface, as it can be evidenced according 
to the model sketched in Fig. ([6li) and invoked in subsection IIV Cl but the calculation 
will not be developed further and the reader interested in this question is invited to refer 
to ref.j^ and references therein. An other residual source of barrier inhomogeneities lies 
in the possible fluctuations of the oxide thickness, because of the link between these two 
quantities as discussed in the next section. Anyway, indications derived from modeling 
the electrical properties of the junctions match very well the conclusions of the HRTEM 
and EELS studies developed in section [TTl which showed the sharpening of the interfacial 
oxide layer after annealing and led to the very probable hypothesis of a chemical mechanism 
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involving bond strengthening between the hrst oxygen layer and the last carbon one. The 
main consequence of this bonds transformation shows out through the very large change 
occurring in the averaged barrier heights of interfaces (a) and (c), with a difference near 
1.4 V, resulting from the partial or complete cancellation of the dipoles 0 ~'''and also 
other dipole changes as discussed in the next section. The case of the diodes set IIB in Fig. 
dH]) does not comply with any previous models, a fact which suggests that the main part of 
the current in these diodes is not governed by an exponential law depending on the applied 
voltage, which was a common hypothesis of all the previous calculations. Therefore, leakage 
currents like those which might occur in threading dislocations or defects appearing as dark 
spots in cathodoluminescence images 28|] could be responsible for a behavior characterized by 
rather higher ideality factors and constant apparent barrier height in Fig. (|H1IIB). Finally, 
homogeneous barrier heights 4*^°™ appearing in Table [TTl can be considered as the most 
representative of the physical situation of each interface, except that unannealed interfaces 
of type (a) may experience two or three different values like mentioned in Table IIIIl the 
lowest being about 0.2 V below that indicated in Table UTl for a homogeneous set of diodes. 
This local barrier lowering may agree the quantity A$ invoked in subsection IIVEI over 
areas with lateral size larger than the Debye length to built the last electrical model relying 
on non pinched-off channels in the previous subsection. It can be concluded that, despite 
some minor imperfect agreements between experimental data and models of inhomogeneous 
barriers detailed in this section and their limitations to the depletion regime of Schottky 
junctions, more informations and deeper physical insight into the situation of these metal- 
diamond interfaces have been gained and can be reinvested in barrier height models discussed 
in the next section. The most important quantity which can be extracted from the analysis 
of electrical characteristics of Schottky diodes is obviously the true potential barrier height 
at an interface free of inhomogeneities and it has been dehned at the beginning of this study 
as It is always underestimated if the barrier height is measured from the slope of a 

Richardson plot and a relevant value needs to be deduced either from adjustment of 
the measured barrier height as a function of temperature, or from the direct adjustment 
of the current-voltage density with the appropriate model chosen among those described in 
this section, or from the statistics of barrier heights as a function of ideality factors in a set 
of diodes, or preferably from several methods. 
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V. BARRIER HEIGHT MODELS AT METAL-OXYGENATED DIAMOND IN¬ 
TERFACES 


The barrier heights issue at metal-diamond interfaces belongs to the long time debated 
topic of how the Fermi level position is determined at semiconductor surfaces and interfaces, 


discussed as early as in the thirties by Sir Nevill Mott 


49| concerning the rectihcation theory 


and by Walter Schottky [50|. Because this first model, which stated the linear dependence of 
the barrier height upon the metal work function with a unity slope, did not comply with the 
huge amount of data gathered since its publication, other explanations have been proposed to 
justify the weaker dependence experimentally observed, known as Fermi level pinning. New 
barrier height models have been progressively introduced within the framework of the ’’metal 
induced gap states” principle [Ml or its generalization to ’’interface induced gap states” [2^ in 
order to link band alignment and either the work function or electronegativity of the metal 
at the broadest possible panel of metal-semiconductor interfaces. All rely on the existence 
of an interface state density localized in the first elementary cells of the semiconductor 
which is neutral when the Fermi level coincides with an energy level <Fo characteristic of 
the semiconductor and referred as to the neutral level, defined here as the energy difference 
between such a level and the valence band edge at interface. Deviations from this position 
and from neutrality occur if a metal having an electronegativity different from that of the 
semiconductor is put in intimate contact with the semiconductor surface, assuming that 
the interface states inside the semiconductor remain invariant. This last hypothesis leads 
to a negative feedback effect due to the reaction of the charge inside interface states onto 
the voltage jump due to the double layer of the areal dipole astride the interface, which 
induces a linear relationship between the barrier height and the difference of the metal and 
semiconductor electronegativities Xm — W, with either a positive slope Sx on a n-type 
semiconductor, or a negative one, —Sx, on a p-type semiconductor as here: 


eo<h|7 = 


X. 


where 


= 


A 


X 


(18a) 


(18b) 


1 + (^0 / ^ o ) -Djs ddi 

with Dis the interface state density and Ax a coefficient which depends on the electronega¬ 


tivity scale, equal to 0.86 if the Miedema scale is used 


24l |. This theory was rather successful 
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for intimate metal-semiconductor contacts, which have been better and better controlled at 
an atomic scale for the seventies, sometimes with an epitaxial metal layer. However, some 
dispersion still remained, which was not attributable to experimental uncertainties. This 
is specially true for diamond, which has not been the subject of extensive reports since 
the nineties, excepted to evidence the wide range of barrier heights existing for the same 
metaljs] and with the notable exception of the hydrogenated surface}]^. Such a dispersion 
has been ascribed to the various structures, terminations, ad-atoms and natures of bond 
which can exist at the diamond surfaces and interfaces, like already invoked in the case 


of oxygen termination. Recently, measurements of barrier heights have been pub 


carefully prepared oxygenated and eventually post-annealed diamond surfaces 


18 


ished at 


52l-l56|. 


Together with the present results for Zr-diamond interfaces, ideal barrier heights are plotted 
in Fig. dH]) with error bars added according to the rules discussed in the previous section: 
either only above the published value when the method underestimated the ideal homoge¬ 
neous barrier or symmetrical when is directly indicated by authors or reassessed 

with the help of an appropriate method among those discussed in the preceding section. 
Despite some dispersion, illustrated sometimes in Fig. (|9]) and (1T0|1 with different barrier 
heights reported for the same metal in literature, two groups of data appear well separated, 
one for the contacts annealed at typical temperatures less than 250 °C or as-deposited and 
an other for contacts annealed at temperatures in the range 400-650 °C, as a function of 
the Miedema electronegativity difference between metal and carbon, the last one being left 
constant at 5.7 eV. Electronegativity was originally used to describe the interaction between 
two (or more) atoms in a molecule and thereafter in solids. But it is impossible here to 
describe the two electronic situations regarding barrier heights in nna nn ealed and annealed 
interfaces with only two electronegativity values, one for the metal and the other for oxy¬ 
genated diamond, since changes after annealing occurred while oxygen still remains at the 
interface as demonstrated in section m Even if the diamond electronegativity were shifted 
by about —1.5 eV to take the oxygen terminations into account in the unannealed interfaces, 
although this is not a very sensible operation, the whole set of data does not show clearly 
any indisputable trend, testifying a hrst shortcoming. A second one appears if the search 
for a linear relationship is restricted to the second group of data, which shows a reasonable 
alignment depicted by the dashed line in Fig. ([9]), but characterized by a positive slope, in 
contradiction with Eq. fllSap and what is experimentally found for all other cases including 
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Electronegativity difference - Xc(eV) 


FIG. 9. Ideal barrier heights as a function of the Miedama electronegativity difference between 
metal and carbon (assumed to be 5.7 eV) for oxygenated interfaces annealed below 250 °C (full 
squares) and for those annealed in the range 400-650 °C (open squares), together with a linear fit 
in dashed line in the second case. The dash-dotted line (H) is the best ht of barrier data on the 

n 

hydrogenated surface by Tsugawa|16l| et al. as a function of the same abscissa. Multiple data for 
the same metal (Au and Mo) come from different publications. 


that of a hydrogenated diamond surface [l6l|. This is the second shortcoming. Incidentally, 
it is difficult to believe that the energy distribution of ’’interface induced interface states” 
may remain the same when going from the hydrogenated then air-exposed diamond surfaces 
to oxygenated ones, as it is well known that the former is provided with adsorbates and 
conducting whereas the later is notj^. An other approach, relying on physico-chemical 
quantities able to influence interface bonding, is necessary. 

In the last twenty years, scientists performed modihcations of the Schottky barrier heights 
by means of various ad-atoms and/or physical, chemical and thermal treatments applied 
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FIG. 10. Ideal barrier heights as a function of the difference between the ionization energy of 
diamond and the work function of metals (multiple data come from different publications) for 
(i) as-deposited interfaces (full squares) on an oxygenated (lOO)-oriented surface, (ii) for contacts 
annealed in the range 400-650 °C (open squares) on the same initial surface, together with a linear 


fit in dottec 
from ref. 


line, and (iii) for contacts on the hydrogenated (lOO)-oriented surface (open circles) 


16l ]. Multiple data for the same metal (Au and Mo) come from different publications. 


The ionization energy includes the electron affinity shift due to surface terminations, yielding (i) 
X -|- Axis) = 1.7 eV; (ii) x + AxiS) = 0.3 eV and (iii) x + AxiS) = —1.1 eV respectively for 


the oxygenated, bare and hydrogenated (lOO)-oriented surfaces from ref. 


331. The dash-dot line (a) 


displays Eq. (jlOap where the last term is ignored, thus featuring the Mott-Schottky relationship. 
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to the surface generally before metal deposition excepted for annealing. They realized that 
these changes are dependent on bonding between interface atoms, and that chemical inter¬ 
actions cannot be ignored between atoms in the neighborhood of the interface. A hrst im¬ 
provement of the theory has been brought by R.T. Tung who showed that potential ener gies 
in the solids and through the interface were more appropriate to describe band alignment 
in such a situation. This idea is applied to data addressed previously which are now plotted 
in Fig. flTOj) . Secondly, the dipole strength of the various double layers occurring at interface 


is the relevant quantity to take the possible interlayer atoms, the different bonds 
even the possible new states due to transfer doping on the hydrogenated surface 


58] and 
16] into 


account. It is lumped into the last term in the following continuity equation which makes 
use of the metal work function (j)M and electron affinity of the semiconductor as follows: 


Co— {Isc{S) — (j)M) — eo dbd Nhd/i£i£o) (19a) 

where the last term is the potential energy jump due to a surface dipole comprising a charge 
per bond Aq on the metal side and —Aq on the diamond side separated by the distance dbd, 
with an areal density of dipoles Nbd, while the ionization energy given by: 


IsciS) = Eg + X + Ax(F) (19b) 

depends on the diamond band gap Eq, electron affinity y of the bare surface and possible 
change in electron affinity Ax{S) due to foreign atoms termination or particular reconstruc¬ 
tion on each surface S, the y -|- Ax{S) values which are retained being indicated in the 
caption of Fig. [101 


When the last term is neglected in Eq. (IIbah , the Mott-Schottky relationship for contacts 
on p-type semiconductors is found, with a unity slope, sketched by a dash-dotted line in 
Fig. (ITOD . If the dispersion of data is forgotten in a hrst approach and the clouds of points 
related to the three categories of interface in this graph are replaced by their respective 


barycenters, one can guess that they align roughly on the Mot 
been already evidenced in silicides on n-type silicon interfaces 


-Schottky line. This trend has 


59( 1 . The inhuence of the cloud 


related to the hydrogenated surfaces should be minored because Ft, An, Cu and Ag metals 
contacts are displayed with zero barrier heights although they should be negative but cannot 
be actually measured and reported. More accurately, if one observes the barrier heights of 
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each metal (Pt, Au, Mo, Al, Zr) whose contact on diamond has been measured in the two 
groups of oxygenated interfaces, respectively before and after annealing, their difference lies 
in the range 1.40 ± 0.25 eV for all of them, corresponding to the same change in electron 
affinity. This common shift indicates that the barrier height for a given metal is primarily 
controlled by the electron affinity of diamond, with a little influence of the dipole due to 
interface bonding, reflected in the last term of Eq. fllOap . and testihes the disappearance 
of the oxygen terminations proper to initial oxygenated surfaces after annealing, or their 
transformation, irrespective of the presence of an oxide or not. As it was discussed in 
subsection IIV Al hydroxyl and carboxyl terminations likely exist before annealing and they 
can participate to bridged hydrogen bonds as suggested in [3^, thus allowing adhesion of 
the atomic layers deposited on the surface by some kind of weak bonding. On the contrary, 
stronger, and therefore shorter bonds between carbon and atoms of the Zirconium oxide 
interlayer, that are mainly oxygen ones, would occur after annealing, in close correlation 
with the shrinkage of the oxide interlayer detected in section [TTl But this statement does not 
mean that ” interface induced gap states” have been suppressed or have never existed. They 
are just accommodating the electric charge which guarantees that the dipole layers, both 
necessary for setting the electron affinity and the bond dipole counterparts on the diamond 
side, can be present. Now, at a more subtle level, the general trend in each group of contacts 
made on the oxygenated surface has an inverted slope in comparison to the Mott-Schottky 
line, as already noticed in Fig. |9l The correlation coefficient is better in annealed interfaces, 
probably because the electron affinity dispersion due to various oxygen terminations has 
been suppressed. The physical reason for the linear relationship depicted by a dotted line 
in Fig. ffTOj) can be thus discussed with more conhdence in the annealed interfaces than in 
unannealed ones because the exact electron affinity is unknown in the second case for each 
metal and therefore replaced by the single value measured in vacuum as indicated in the 
caption of the hgure. First, it is useful to notice that neither the hrst term in Eq. (I19ap . which 
also scales the abscissa in this diagram, nor the second term can account for the observed 
trend if this later term is supposed to be proportional to the former one because the slope 
would be increased above one when the ionization energy is held constant. Therefore, one is 
forced to conclude that hidden physical parameters able to govern this second term might be 
involved. These are the nature of the bonds and interlayer, the charge transfer across them, 
their polarizability, and the interface structure with a special attention to the interlayer 


48 









thickness and relative permittivity. In the case of metals like Pt, Ir and An, oxide formation 
is very unlikely since the enthalpy of formation is either positive or close to zero, whereas 
an oxide is present when Zr, as shown at the beginning of this study, and most probably 
Al, are used. If the interlayer is supposed to be a dielectric consisting of a ionic solid, its 
thickness tox is added to the bond length du so that the total dipolar distance is roughly 
dbd = tox + du. Then this distance determines the strength of the last term in Eq. fll9ap 
since the charge counterpart of the dipole on the metal side is now located at a slightly 
larger distance than tox from the last diamond atomic plane and if the charge inside the 
oxide layer is neglected at first approximation. Moreover, important effects are due to the 
charge transfer Aq which is governed by the nature of the bonds, electronegativity difference 
between the metal and semiconductor, the polarizability of bonds and interlayer thickness, 
as described in the bond polarization model proposed by R.T. Tung and also screening 


which can mainly accounted for by the relative permittivity occurring at the denominator 
of the last term in Eq. (I19ap . To get a deeper insight in the problem at least at a first 
order level of approximation, let us define the abscissa Wsm = {Isc{S) — 0m) and take 
the derivative of the barrier height cq when the electron affinity, and hence ionization 
energy are constant: 


d{eo 


= 1 + 


d{eo Aq duNu/iei Eq)) 


( 20 ) 


dWsM d(j)M 

This slope would stay close to one whenever the second term is negligible but can become 
negative as observed in Fig. ffTOj) if both the second term is negative and its absolute value 
exceeds unity. Since the main trend occurring when the metal is changed to another one 
with an increased work function involves a decreased reactivity with oxygen and therefore 
a thickness reduction or disappearance of the oxide layer, one must wonder what are the 
resulting effects on the second term of Eq. fl20l) . even if a continuous variation is more artificial 
than reality. The dipolar distance d^d is obviously decreased while the interlayer permittivity 
Ei initially near that of high-k dielectric materials is going down closer to unity, resulting 
in a smaller screening effect. But because theses two parameters are following the same 
direction of variation and influence the second term by their ratio, the resulting sign of the 
derivative turns out to be uncertain. Consequently, the fundamental issue relies mainly on 
the charge transfer Aq variation, which is dependent on the effective bonding at the metal- 


semiconductor interfaces 


2d, 


3l| . including not only orbitals of the two original materials 
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but also atoms or molecules belonging to the interlayer. In such a molecular orbital scheme, 
metals with the lowest work functions (and also more electropositive) would make more ionic 
bonds, with a polarizability likely augmented by the presence of additional oxygen bonds 
inside the interlayer, thus inducing a stronger positive charge transfer Aq in the orbitals 
of the last metal atomic layer. Conversely, noble metals (and also more electronegative 
ones) are expected to share the electrons of their d states in a more equalized way with 
carbon orbitals, while carbide forming metals would follow a medium behavior. Such a 
general trend of the charge transfer Aq which would experience a decrease when the metal 
work function (and electronegativity) is increased, is able to justify the systematic negative 
sign of the derivative in the second member of Eq. fl2Up . Consequently, the slope of barrier 
heights as a function of Wsm = — 4>m) can take values smaller than one, as 

often observed in many semiconductors, and in the case of a magnitude of the absolute 
value of the last term in the same equation greater than unity, it can become negative as 
shown in Fig. (ITOll . in contrast to what is known in all other semiconductors. Such a model 
can rely on the electronegativity difference between the semiconductor and metal, which 
mainly governs the dipole strength induced by interface bonding, but only for the last term 
in Eq. fllQap . independent of the interface induced gap states at hrst order. At a second 
order level of approximation, the bonding nature of atoms inside the oxide and fixed charges 
in it will be also able to bring changes in Aq and these quantities probably contribute to 
the small deviations visible in Fig. flTOl) for annealed interfaces, together with the oxide 
interlayer thickness, permittivity, and changes in the atoms arrangement at interface and in 
the hrst diamond atomic layers. The same model also applies both to unannealed oxygenated 
interfaces and hydrogenated interfaces, however with larger huctuations in comparison to a 
linear behavior, and a positive slope in the second case. It covers therefore all the situations, 
although it is difficult to afford quantitative parameters in the second term of Eq. fll9ap for 
each metal in an accurate way. 


VI. CONCLUSION 

In this work, interfaces of Zirconium with oxygenated diamond surfaces before and af¬ 
ter two kinds of anneal has been characterized on the one hand by HRTEM imaging and 
EELS spectroscopy, and on the other hand by thorough electrical measurements made on 
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the Schottky junctions. HRTEM and EELS techniques evidenced an oxide interlayer com¬ 
prising approximately one or two Zr oxide atomic layers, which sharpened after annealing. 
Electrical measurements were performed to collect and analyze current-voltage character¬ 
istics at several temperatures in many diodes in order to investigate the detailed shape of 
the current density behavior with applied voltage, the Richardson constant, the apparent 
barrier height and ideality factor. Models taking the barrier height and current density 
inhomogeneities into account have been implemented and brought face to face with experi¬ 
mental results. From these comparisons, several physical parameters have been estimated. 
The Richardson constant turned out to be close to one tenth of the standard value, and 
even one thirteenth when the barrier height disorder reached its maximum value. The ef¬ 
fect of different shapes of barrier fluctuations has been studied either numerically in the 
case of local barrier enhancements or from two analytical models adapted to local barrier 
lowering and current crowding in the active layer of the diode. In the former case, which 
happens when fixed repulsive charges are close to the interface like evidenced in metal-Zr02 
-diamond structures, the image force effect is the only relevant one as far as current-voltage 
characteristics are concerned. In the case of random patches comprising a smaller barrier 
than the ideal homogeneous one, it has been shown that barrier inhomogeneities can be 
characterized by several parameters related to their amplitude, density, and to the pres¬ 
ence or absence of both the pinch-off effect in these lowered barrier patches and current 
crowding in narrowed channels of the depletion zone. Mixing two of the last effects turns 
out to be necessary to quantitatively explain the curvature of current-volt age characteristics 
of unannealed interfaces in a semi-logarithmic plot. The ideal homogeneous barrier is the 
truly representative physical quantity of the interface because it is free from the local and 
random inhomogeneities which are responsible for the non ideality of the diodes. It has 
been inferred from matching previous models to electrical data. Rules and uncertainties 
have been derived for a confident evaluation of this fundamental quantity. All other param¬ 
eters have been extracted in turn from the comparison of the non ideal behaviors with the 
quantitative predictions of the models. Main conclusions were the strong reduction of both 
the barrier inhomogeneities and homogeneous barrier of 1.4 eV after annealing interfaces at 
450 °C. Because of the known inhomogeneity of the oxygen terminations and enhancement 
of the electron affinity by also 1.4 eV on the initial oxygenated diamond surface, these facts 
strongly suggested that original oxygen terminations vanished after such an annealing and 
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are replaced by stronger bonds between carbon and atoms of the Zirconium oxide interlayer. 
Then, the more general issue of barrier height determination and Fermi level pinning has 
been addressed by means of an analysis of the experimental data collected in the present 
work and in recent literature, reassessed with the help of the methods presented before 
when necessary. The same evolution of barrier heights for four other metals after anneal¬ 
ing conhrms the disappearance of oxygen terminations, irrespective of the presence of an 
oxide interlayer, and the electron affinity of the diamond surfaces as a relevant quantity for 
barrier height description. For the whole set of metals in contact with either unannealed 
or annealed oxygenated surfaces, or hydrogenated surfaces, a rough but clear alignment of 
barrier heights along the Mott-Schottky straight line, with a unity slope, really emerged 
regarding the difference between the ionization energy of diamond surfaces and metal work 
function. However, for each particular interface, the behavior at a smaller scale appeared 
different. In the case of annealed oxygenated interfaces, the slope turns out to be slightly 
negative, a result incompatible with the interface induced gap states model which assumes 
these states as invariant from one metal to another. On the contrary, the bond-polarization 
scheme, which relies on the interface dipole strength, itself dependent of the molecular or¬ 
bitals involved in the bonds across the interface, may justify the slopes both at a global level 
and for a particular surface. All these concerns have been presented and discussed with 
the help of new experimental results, improved models and recent high quality experimental 
works in the purpose of fostering new advances in understanding the old problem of potential 
barrier at metal-covalent semiconductor interfaces. 
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